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This paper describes the use of an Eulerian Dispersed Phase (EDP) model to simulate 
the water injected from the flame deflector and its interaction with supersonic rocket 
exhaust from a proposed Space Launch System (SLS) vehicle. The Eulerian formulation, as 
part of the multi-phase framework, is described. The simulations show that water cooling 
is only effective over the region under the liquid engines. Likewise, the water injection 
provides only minor effects over the surface area under the solid engines. 
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constant pressure specific heat of fluid, J/kg-K 
constant pressure specific heat of particle species z, J/kg-K 
average particle diameter of particle species z, m 

total internal energy (thermal and kinetic energy) of particle species z, J /kg 

Nusselt equation correction factor of particle species z 

buoyancy force on particle species z, N 

interphase drag force on particle species z, N 

lift force on particle species z, N 

pressure gradient force on particle species z, N 

turbulent dispersion force on particle species z, N 

virtual mass force on particle species z, N 

Latent heat of evaporation of species z, J/kg 

rate of production of mass for particle species z, kg/m 3 -s 

number density of particle species z, m " 3 

pressure of particle species z, Pa 

Prandtl number 

saturation pressure , Pa 

reference saturation pressure, Pa 

ambient temperature, K 

fluid temperature, K 

temperature of particle species z, K 

evaporation temperature of particle species z, K 

vector velocity of particle species z, m/s 

heat flux per particle species z, W/m 3 

gas constant of the evaporation of particle species z, m 2 /s 2 -K 
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ri average particle radius of species z, m 

a v correction factor in the Hertz-Knudsen equation 

\f heat conductivity of fluid, W/m-K 

[if dynamic viscosity of fluid, Pa-s 

e P i emissivity coefficient of particle species z 

a Boltzmann constant of radiation, 5.67 x 10" 8 , W/m 2 -K 4 

p P i “fluidized” density of particle species z, kg/m 3 

Ppi material density of particle species z, kg/m 3 

Tn temperature relaxation time of particle species z, s 

T U i Stokes relaxation time of particle species z, s 


I. Introduction 

Since the retirement of the Space Shuttle, NASA has been developing a new heavy-lift capability, the 
Space Launch System (SLS). In conjunction with this development, NASA is re-assessing their launch site 
capabilities. One key component of particular interest is the Main Flame Deflector (MFD). The main 
purpose of the MFD is to minimize the plume impingement effects on the launch facility and to minimize the 
induced environmental effects on the vehicle. The MFD must be protected from harsh thermal environments 
generated by the plume impingement. 

Currently, the MFD is kept cool by water sprays from the nozzles mounted on the apex of the deflector 
(Figures 1). The MFD water is one of the many water suppression systems employed at the launch complex. 
There are also water nozzles mounted around the exhaust hole to reduce the ignition overpressures and water 
rainbirds mounted on the launch deck to suppress the launch acoustics. This study focuses only on the MFD 
water nozzles, which can deliver water flow rates up to 300,000 gallons per minute (gpm) on both sides of 
the deflectors (Figure 2). 

Recent advances in computational fluid dynamics (CFD) technique and computing resources have pro- 
vided the necessary predictive capabilities during the MFD redesign. 1,2 However, these predictions often 
neglect multi-phase modeling to avoid model complexity and to reserve conservatism associated with un- 
cooled or dry rocket plumes. Since the conservative results show severe thermal environments, it is prudent 
to know the effect of the water suppression system on the MFD in order to ensure that the MFD is not 
overdesigned. 

The interaction of water sprays with supersonic flow is a complex phenomenon as it involves condensation, 
evaporation, diffusion, and dispersion processes. In this paper, the Eulerian Dispersed Phase (EDP) model 
is employed to model the water sprays into the gaseous environment. The EDP is coupled with the multi- 
species gaseous flow via source terms that govern the mass, momentum and thermal transfers between the 
two phases. The solutions allow designers to evaluate the effectiveness of water suppression systems on the 
MFD. 


II. Modeling 

The analysis pertains to a complex system involving multi-phase, multi-species, and complex geometry. 
These complexities provide a significant computational challenge. All modeling reported in this paper was 
performed using CFD+4- 12.1, a multi-purpose CFD software suite developed by Metacomp Technologies, 
Inc. 3-5 

A multi-block, unstructured mesh is created to model the two solid and four liquid engine configurations 
on the launch pad (Figure 3). The mesh comprises of 34.3 million cells. A symmetry boundary condition 
is used to reduce the computational domain by half. The entire domain is initialized with a quiescent flow 
field. 

The solution is based on compressible, multi-species, Reynolds-Averaged Navier-Stokes (RANS) with a 
realizable k-e turbulence model. The multi-species simulation utilizes four species groups: SSME, SRB, 
water vapor, and air. The properties required to describe these groups include gas constants, molecular 
weights and specific heats. For the SSME and SRB groups, those properties were obtained from running a 
chemical equilibrium with applications (CEA). 6 The simulation starts from the combustion chamber where 
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the stagnation pressure and temperature are specified ahead of the throats. Standard non-reflecting, inter- 
polated boundary conditions are imposed at the farfield boundaries, as well as no-slip conditions at solid 
surfaces and a symmetry boundary condition at the midplane. Steady-state solutions are obtained using an 
algebraic multi-grid accelerated relaxation scheme. 

The crest water injected from the apex of the MFD is modeled as a rectangular block with 8 slits along 
one side of the block for the water inflow boundaries (Figure 4). The flow areas, inflow velocity, and the 
canting of the nozzles are all modeled for a better representation of the collective jet. The actual water 
nozzles are canted 5 degrees above and below the horizontal axis, alternately. These angles of flow direction 
are specified in the water inflow boundary conditions. 

III. Governing Equations 

The general form of conservation of mass, momentum, and energy equations that describes the particles 
in the dispersed phase and the fluid in the continuous phase in a coupled manner as follows: 

■ + v • (Ppiiipi) = rhi 

+ V • ( PpiUpiUpi ) = F Dl + FvMi + Ftdi + Fl, + FpGi + F Bt + m l u pi 
9(ppi e pi) ( f ^ o \ _ A I 

Qt 4“ V • (PpiUpiUpi) — Qpi "b ^i^pi 

where Fp is the interphase drag force, Fvm is the virtual mass force, Ftd is the turbulent dispersion 
force, Fl is the lift force, Fpc is pressure gradient force and Fp is the buoyancy force. 

The interphase drag force is given by: 


where the Stokes relaxation is: 
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The heat transfer term is given by: 

Q pi = P V iC Vi ^ (Tf - T pi ) + 6 £ f P fa (T a 3 4 m6 - T 4 ) 

T Ti Mi Ppi 

The first term arises from the heat transfer between the fluid and the dispersed phase. The second term 
is the radiative heat loss of the particles where T am b is taken to be zero. The thermal relaxation time is 
defined by: 


where, 




Pr _ aj£pf 

x f 

The evaporation source term is written as: 

rhi = Nippier i) 

The evaporation rate is solved by the Boiling and Hertz-Knudsen models. The evaporation rate for the 
Boiling Model is given by: 
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The evaporation rate is solved by the Hertz-Knudsen model for evaporation of volatile liquids by: 
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IV. Discussion of Results 

In order to capture the physics associated with water jet injected into the supersonic plumes, 34.3 million 
cells with mesh refinements in the jet-plume interaction region and prism layers on the trench surface are 
constructed (Figures 5-6). Figure 7 shows the Mach number contours at the planes cutting across the SRB 
and SSME nozzles. It can be seen through the plume exhaust that the SRB nozzle is slightly under-expanded 
while the SSME nozzles are clearly over-expanded. Near the top of the deflector, the plume is slightly altered 
due to water injection and the same observation can be made for the temperature contours in the same plane 
cuts, shown in Figure 8, where temperature is in unit of Kelvin. It is expected that little to no water is 
present at locations where plume temperature exceeds 373K; this phenomenon is confirmed in Figure 9, 
which shows the water vapor species at the SRB and SSME cutting planes. The temperature contour on 
the MFD looking from the top is shown in Figure 10. It should be noted that the water injection angles are 
alternately canted 5 degrees up and down from the horizontal axis. The fact that these patterns are not seen 
on the deflector surface indicates that water quickly vaporizes before it reaches the surface. It is apparent 
that the cooling effect is much greater in the region under the SSME and water injection only reduces surface 
temperatures on the other side of the MFD, under the SRB. This behavior can be explained by observing 
the EDP density at a plane cut near the SRB (Figure 11), which reveals that the SRB plume pushes the 
water jets over to the other side of the MFD. 

Furthermore, a comparison of surface temperature between the simulations with and without water 
injection is shown in Figure 12. The temperature is plotted along the surface of the MFD at a plane 
cut across the SRB. For the coarse- mesh simulations the temperature drops drastically near the water jet 
nozzle and also along the other side of the deflector, where there is little to no presence of the plume. 
Likewise, water has little to no cooling effect over the MFD surface where the plume is deflected along. The 
explanation to this is because the SRB plume temperature exceeds the water evaporation temperature by 
several orders of magnitude, so when water is injected to the plume; it quickly vaporizes before reaching the 
surface of the MFD. The water-on simulation with fine mesh and prism layers shows some improvement on 
the deflector surface temperature, but it is still much higher than water evaporation temperature over these 
plume impingement areas. The same temperature drop is observed near the water nozzles for coarse and 
fine mesh simulations with water injection. 


V. Conclusion 

A multiphase simulation was carried out to determine water cooling effects on supersonic plume impinge- 
ment, in an effort to study the plume-induced environment due to SLS launch. It has been shown that 
cooling effect greatly depends on how the vehicle is positioned over the MFD. If water is injected directly 
onto the rocket plume, it would evaporate quickly and has little effect on the MFD. The water cooling would 
have been more effective if it were installed along the MFD, as shown in another study. 7 

It is possible that the water jet angle can also have a cooling effect. For that reason, the water nozzles 
might need to be canted to a larger angle downward. This might pose a challenge in the MFD design, but 
it is an option to be considered. 
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Figure 1: Launch Complex-39 Main Flame Deflector. 


Figure 2: MFD Water Nozzles. 
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Figure 3: Truncated SLS vehicle above the Main Flame Deflector. 



Figure 4: MFD Apex. 
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Figure 5: Grid refinements in the flame trench. 



9 of 13 


American Institute of Aeronautics and Astronautics 




M 



M 

S 4.84470400 
3.87570400 
2.90686400 
1.93790400 
9 . 68930-01 
0.00000400 




Figure 7: Mach number contours in the SRB and SSME cross sections. 
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Figure 8: Temperature contours (unit of Kelvin) in the SRB and SSME cross sections. 
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Figure 9: Water vapor mass fraction in the SRB and SSME cross sections. 
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Figure 10: Surface temperature contours (unit of Kelvin) on the MFD (top view). 
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Figure 11: EDP density near the water jet and SRB. 
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Figure 12: Surface temperature near the SRB cross section. 
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